Implementation of many-qubit Grover search with trapped ultracold ions 



00 

o 
o 

Of 

D ' PACS numbers: 

oo 



Oh 



oo 
o 



X 



Wan-Li Yang 1,2 , Hua Wei 1 ' 2 , Chang- Yong Chen 3 , and Mang Fengfl 

1 State Key Laboratory of Magnetic Resonance and Atomic and Molecular Physics, 
Wuhan Institute of Physics and Mathematics, Chinese Academy of Sciences, Wuhan 430071, China 

2 Graduate School of the Chinese Academy of Sciences, Bejing 100049, China and 

3 Department of Physics and Information Engineering, 
Hunan Institute of Humanities, Science and Technology, Loudi 417000, China 

We propose a potentially practical scheme for realization of an n-qubit (n>2) conditional phase flip 
(CPF) gate and implementation of Grover search algorithm in the ion-trap system. We demonstrate 
both analytically and numerically that, our scheme could be achieved efficiently to find a marked 
state with high fidelity and high success probability. We also show the merits of the proposal that 
the increase of the ion number can improve the fidelity and the success probability of the CPF 
gate. The required operations for Grover search are very close to the capabilities of current ion-trap 
techniques. 
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I. INTRODUCTION 



Grover's quantum search algorithm [1] has been considered to be an efficient amplitude-amplification process for 
quantum states by exploiting the parallelism of quantum mechanics. As a remarkable idea in quantum computation, 
Grover search algorithm could effectively exemplifies the potential speed-up offered by quantum computers. Recently, 
many proposals have focused on the search algorithm by the adiabatic evolution method [2-4] or by nonadiabatic 
scenario [5] . Although achievement of these schemes needs stringent conditions and demanding techniques [6] , they are 
really wonderful ideas. On the other hand, some authors had addressed the effect of decoherence [7], gate imperfection 
or errors [8], and noise [9] on the efficiency of quantum algorithms. We also noticed that, there had been intensive 
■ interests in achieving Grover search algorithm theoretically and experimentally by using NMR systems [10], linear 
CT) optical elements [11,12], trapped ions [13-15], cavity quantum electrodynamics (QED) [16-19], and superconducting 
. mesocircuits [20]. 

Trapped ions have been considered to be a promising candidate for quantum-information processing (QIP), due to 
long coherence time of qubits, full controllability of operations and high efficiency of detection. In [13], a many-qubit 
0^ ' Grover search algorithm based on the hot-ion quantum computing [21] in decoherence-free subspace was proposed, 
where collective dephasing errors could be kept away from the qubit encoding. Ref. [14] is an alternative for simulating 
Grover search algorithm in an ion-trap system, in which the search could be carried out by more than two qubits. 
However, in order to achieve the conditional phase flip (CPF) gate, each ion should be illuminated by two lasers with 
different polarizations. In addition, the proposal needs not only elaborately designed sequences of lasers, but also 
the auxiliary states of the ions, which would result in considerable difficulty in experimental realization with growing 
number of the qubits. 

In this paper, we propose a potentially practical scheme for implementing n-qubit (n > 3) Grover search in an ion 
trap. Compared to above mentioned schemes, our scheme includes following merits: We carry out the CPF gates 
with nearly unity success probability and fidelity in a straightforward way, which could relax the rigid requirement 
on accurate sequences of laser manipulations in previous work, such as sequences for one-qubit and two-qubit gates 
[12,22] or swap gates [23]. Moreover, the increase of qubits in our scheme could improve the fidelity and the success 
probability of the CPF gate, which are favorable for a scalable Grover search. 

Section II describles the general method for a CPF gate in a linear ion trap. We then describle, in Sec. IV, the 
implementation of n-qubit Grover search algorithm based on the proposed CPF gates, which is almost within the 
reach of the present technology and extendable to other QIP candidate systems. We conclude with a discussion in 
Sec. V. 
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II. N-QUBIT CPF GATE 



We consider n identical three-level ions confined in a linear ion trap, as shown in Fig. 1 where the ions are 
individually irradiated with traveling wave laser fields tuned to the first lower vibrational (i.e., red) sideband. The 
three-level ionic states under our consideration are denoted by \ fj) , \gj) ,and \ej} , with \ fj) and \gj) being states lower 
than \ej). Because the resonant transition happens between \gj) and \ej) by the laser field, \fj) is not involved in the 
interaction throughout our scheme. In the rotating wave approximation, the Hamiltonian in units of ft = 1 reads, 

n 

H = ^{i/a+a + woo^- + [\E+{r,t)a+ +H.c.}}, (I) 

where E + (r, t) — E a cxp[— i(kz — curf + (f>)] is the positive frequency part of the driving laser with amplitude E and 
frequency lui = loq — v. a + (a) is the creation (annihilation) operator of the center-of-mass vibrational mode commonly 
owned by all the ions, er+ = \ej) (gj\ , aj = \gj) (ej\ , and a z j — ^ C I e j ) ( e j\ ~ are the raising, lowering and 

inversion operators for the jth trapped ion, respectively, v, u>o and A are the trap frequency, the atomic transition 
frquency and dipole matrix element (assumed to be real as convention), respectively. In the resolved sideband limit 
and in the interaction picture, the Hamiltonian of Eq. (1) excluding the terms associated with carrier transition can 
be simplified to [24], 

^/ = EE n i W e " r ' 2/2 ^<i^TiTT (a+) ^ +1 + Hc -> (2) 



where r\ — k/\/2nvM is the Lamb-Dicke parameter with M being the mass of the ion. We assume the laser pulses to 
be time-dependent in Gaussian type, and the time varying Rabi frequency flj(t) of the laser field is given by 

= nL«exp{-(i - t ) 2 /2r*}, (3) 

with Tj being the duration of the Gaussian shaped pulse irradiating on the jth ion. In the Lamb-Dicke regime, the 
Hamiltonian can be approximated by the expansion to the first order in 77, 

n 

H 1= ^2 irfiljWa+aje-*' - aa+e*'). (4) 
i=i 

Furthermore, by assuming <f>j = w/2, we reduce Eq. (4) to 

n 

H 2 =Y^ Vttj (*) {a + °j + oaf). (5) 
i=i 

We first assume that the center-of-mass mode of the ions is initially in the vacuum state |0) , and except the last ion 
(i.e., the nth ion) initially prepared in the excited state |e„) and the fcth ion in the state \gu) , other ions are initially 
prepared in the state |/). So we only have the nth and the fcth ions interacting with the vibrational mode by following 
evolution, 

n-1 



nj =1 1 , j¥fe \9k) \fj) \e n ) |0> — exp(-i J H 2 dt) J] \9k) \fj) K) |0> 



j-J^k 

= [§2 cos(^L) + ^jfi] x IJ \g k ) \fi) \e n ) |0) 

k k j-j^k 

a a ™— 1 

+ ^ L [cos(^)-l]|e fc )| 5n ) J] l/i)l°> 

k 3,3=£k 

-^sinM fe )|. 9fe )|.g„) J] (6) 

k j,j¥=k 
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where ti> k = V^T+ and 



dj= ^ max exp{-(t-to) 2 /r 3 2 K 
Jo 

= [° 04 ax exp(-t 2 /2r2)rfi 

■/-to 



«^ max V2^-, (7) 

where erf [z] = (2/ v / 7r) J q 2 e _t dt is the error function, and j = 1, 2, • • • , n. In concrete calculations, we assume that 
erf [to/v/^Tj] — ► 1. Next, we consider another situation, that is, the last ion is initially in the excited state \e n ), and 
for other ions some of which are initially in the state \g) and the rest are in the state |/), then we can acquire the 
corresponding time evolution, 

11",=! \fi) K) |0) exp(-i / H 2 dt) ]\ \ 9j ) \,n) \e n ) |0) 

= & CO8{r ) + fj=A ]x JJ \ gj) \ fl) \ Bn) |0) 

q s n— 1 

+ ^[cos(^)-l]2<? fc |e fc )| 5n ) l*)l//)|0) 
fe=i j,i=i,j& 

hi ■ ^ 

__l sinM) JJ b)|/ ; )|.g„)|l), (8) 



where $ = yv 2 + Sj=i w i tn s t^ 10 num bcr of the ions initially in the state \g) and other denotations are defined 
as the same as above. To achieve our aim, we need to consider another case: If the ions are initially in the state 
rij=i \fj) l e «) |0) > then the ions, except the last one, do not interact with the vibrational mode. The evolution of the 
system is as follows, 

= [cosM„) \e n ) |0> -isinfa0 n ) |.g„) |1)] J]. =1 l/i) • (9) 



Based on Eqs. (6)— (9), we can construct an n-qubit CPF gate. In our proposal, the qubit definitions are 
the same for other ions, except the last ion, i.e., the logic state |1) (|0)) of the ith qubit is denoted by (\gi)) 
of the ith ion with i = 1,2, ••• ,n — 1, whereas the logic state |1) (|0)) of the nth qubit is represented by \e n ) 
(\g n )) of the nth ion. By considering the quantum information encoded in the subspacc spanned by the states 
i\gi) , l/i) j Iff2> , I/2) , • • • , \g n -i) , |/n-i) , Iff™) , \e n }}, we have Eq. (9) in the case of r]$ n = ir to be, 



IIli M K) |o) — - IT'I \fj) K) |o) . (10) 



r n— 1 t — r n— 1 

[ J=1 i/i>ie->io>— -n i= 

Furthermore, we assume that coupling strenghths satisfy the following condition, 

<?! =0 2 = ...=<?„_! »0 n . (11) 

After inserting Eq. (7) into Eq. (11), we have, 

m =£ = ^<l, (12) 

with z = 1, 2, • • • , n — 1, which implies that the condition in Eq. (11) could be met by adjusting the pulse widths 
Ti and the maximum coupling strength Sl^ax- To keep the ions in the vacuum center-of-mass mode, we assume that 
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each laser pulse width has the identical value (jj = To), but we have different fi^ax by setting m = ^ ax /ff max <C f 
to meet the requirements. Then Eqs. (6) and (8) can be reduced to 



n-l 



Yljl^WMMM^P* II \9k) \fj) \e n ) |0) 



and 



n-l 



U ; / ^ Ift) I//) K) 10) — a s x J] \ 9j ) |/,) |e„) |0) 



(13) 



(14) 



where a s = [to 2 cos(V?ti 2 + sn/m) + s]/(m 2 + s) and /3 = [to 2 cos(Vni 2 + l7r/m) + 1]/(to 2 + 1) = a\. So if the 
requirement of Eq. (12) is met, we have (3 ~ a s ~ 1 and thereby obtain an n-qubit CPF gate in our computational 
subspace, where the ions from the first to the (n — l)th represent the control qubits, and the last ion represents the 
target qubit. We take n — 4 as an example below. From Eqs. (10)-(14), an approximate four-qubit CPF gate can be 
reached as follows: 



U, 



(4) 



CPF 



diag{\, 1, a 3 , a 2 , 1, 1, a 2 , /3, 1, 1, a 2 , (3, 1, 1, /3, -1}, 



(15) 



in the computational subspace spanned by {|5i525354) , \g1g2g3e4} , I5152/354) , \g1g2f3e4) , \g1f2g3g4} , \g1f2g3e4} , 
I51/2/354) , \g1f2f3e4) , I/1525354) , |/i5253e 4 ) , I/152/354) , 1/152/364) , I/1/25354) , |/i/253e4) , 1/1/2/354) , 1/1/2/364)}, 
where a 2 = 0.99952, a 3 = 0.99516, (3 = 0.99988 in the case of to = 0.1. 

We now turn to calculation of the fidelity and the success probability according to the relations F = 

(*o| U ( cl% [*/) (*/| U^l F |*o) [25] and P = {¥} \V f ) , where |*/) is the final state after the n-qubit CPF gat- 
ing has been made. The overline indicates average over all possible components in |^ ) • I n the case of n qubits, 
we set |*o) = -7= (!.-?!> + |/i))(|52) + I/2)) ■ ■ ■ (\g n -i) + |/„_i))(|5„) + \e n )), the infidelity and the success probability 



which gives a general assessment for our gating are obtained straightforwardly by 



Infiderity = 1 - F = 1 



E:= 2 1 d n) a s + (n-l)/? + 2"- 1 + l] 2 
2 n E"= 2 Ci n) a 2 + (n - l)/? 2 + 2™-i + 1] : 



and 



P = 



where we denote the number of a s (s = 2,3,- 
the equations, 



2" 

• , n — 1) by C'i™' 1 with n the number of the qubits, and C, 
U s + ° s +i — °s+l > 



(16) 
(17) 

{n) fulfilling 



and 



S"-2 C ^ = T 



n. 



(18) 



(19) 



We have listed some examples in Table I for the values of C, 



(n) 



TABLE I. List of the values of Q 



(n) 



qubit number n 


number of a 2 


number of a 3 


number of 


number of as 


number of ae 




number of a n -i 


Y^™-1 s-i(n) 
2^s=2 °s 


n = 4 


3 


1 

















4 


n = 5 


6 


4 


1 














11 


n = 6 


10 


10 


5 


1 











26 


n = 7 


15 


20 


15 


6 


1 








57 



We have plotted the infidelity and the success probability versus m — ^max/^max m Figs. 2 and 3, which clearly 
indicate that our proposed n-qubit CPF gate is of high fidelity and high success probability as long as the value of to 
is small enough. The figures also show the increase of F and P with the number of the ions. To obtain the maximal 
fidelity and success probability, provided that Sl l max = Q m , we have made some numerical calculations, which show the 
values of m to be to = ri^ ax /fi m = 0.0122, 0.0122 and 0.0147 corresponding to n —3, 6 and 9, respectively. One can 
see from the figures that the implementation of the n-qubit CPF gate with high fidelity and high success probability 
can be realized by suitable laser-ion coupling strength ratio of the last ion to the others. 
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III. N-QUBIT GROVER SEARCH 



In this section, we will implement an n-qubit Grover algorithm by our gating discussed above. One of the key steps 
in Grover search algorithm is to find a marked element in an unsorted database of size N. Generally speaking, the 
Grover search algorithm consists of three kinds of operations [13]. The first one is to prepare a superposition state 
I^q) = J2f=o 1 K) using Hadamard gates H® n (n being the number of qubits). The second is for an iteration 

including following two operations: (a) Inverting the amplitude of the marked state \p) using an operator 
J p = I — 2\p) (p\ , with / the identity matrix; (b) Inversion about the average of the amplitudes of all states using the 

diffusion transform D^ n \ with = j| — 8ij = 1,2,3, •• • , N) and N = 2". This step should be carried out by 

at least tx\J~N /4 times to maximize the probability for finding the marked state. Finally, a measurement of the whole 
system is made to get the marked state. 

As defined above, the logic state |1) (|0)) of the ith qubit is denoted by |/j) (\gi)) of the ith ion, with i = 1,2,--- , n—1, 
whereas the logic state |1) (|0)) of the nth qubit is represented by \e n ) (\g n )) of the nth ion. By considering the 
quantum information encoded in the subspace spanned by {\gi} , , \g 2 ) , I/2) , ■ ■ ■ , \g n -i) , \ fn-i) , Iff™) , |e„>}, we 
have following transformation, 

g(n) = W® n J$.. W® n J p = W^ n J^l 2 ... gn W® n J p = -D (n) J p , (20) 

where Jqo?..o = 

diag{-l, 1, • • • , 1} = - 2 |00 • • • 0) (00 • • • 0| (the number of is n) and the Hadamard gate in our 
whole computational subspace is given by 



n , 1 

i—1 ^ v 



1 1 




1 1 




1 1 


1 -1 


® 


1 -1 


<s> - - * <s> 


1 -1 



(«) 



In the case of n qubits, Eq. (20) implies that the diffusion transform iM") = — W® n J, 
no matter which state is to be searched. The only change is the operator J p for different marked states. Based on 



(21) 



, W® n is always unchanged 



the CPF gate constructed in last section, we have 
the gate Jqq1. and other relevant operations, 



diag{l, 1, 



<^00--0 ~ tJx.nSx.n-l ' ' ' S x ,lj\\!..\S x ,l 1 1 1 Sx,n— l°x,t 



-1} 



Uqp F , from which we could achieve 



and 



f)(n) = lf®»jW ^«» 



(22) 



(23) 



where S x ,j = \fj)(gj\ + \9j)(fj\ with j ^ n, and a x . n = \e n )(g n \ + \g n )(e n \. To achieve Q^ n \ we will construct the 
CPF gate J p = I — 2 \p) (p\ . For example, in the case of n = 4, the number of possible quantum states is 2 4 , and the 
operation is to label a marked state by J p with p one of the states from {|0000) , 1 0001) , |0010) , • • • , 1 1111)}. To carry 



out the four-qubit Grover search, we need two four-qubit Hadamard gates W® 4 . Based on the gate Jim 
(See Eq. (15)) marking the state 1/1/2/364) , we could construct other fifteen gates for the marking job as, 



U, 



(4) 



CPF 



J 9192 93 34 


= Joaoa 


= ^x 


4&X 


^S x ^S x sJii\iS x ^iS x ^S x ^a x ^, Jg 1 g 2 g 3 e i — J0001 — S X ^S X ^S X ^ J\u\S x ^iS x flS x ^, 


^ '91 92 f 3 94 


= Joaw 


= cr x 




2S x sJl\l\S x sS x ^O x ,<i, Jg 1 g2f 3 e i = 


: JoOll = S Xi 2S Xi i JuuS x ,lS x ,2, 


J 91 f 2 939 A 


— Jowa 


= cr x 




3<Jx,l JllllS Xy iS Xy 3<J x ^, Jgi/293e 4 : 


= JoiOl — S x ,3Sx,lJllllSx,lS x ,3i 


Jgihhgi 


— Joiw 


= cr x 


iS x 


1 JllllO Xj lCr Xi 4, Jg 1 f 2 f 3 e 4 — Joill = 


z S Xi iJimS Xi i 


Jfl 92 93 34 


— J1000 


= cr x 


iS x 


3S x ,2JllllS Xt 2S Xt 30 Xt 4, Jf 1 g 2 g 3 e i = 


JlOOl = <Sx,3<Sx,2 JllllS Xl 2S x ,3, 


Jfl 92 f 3 94. 


= J1010 


= Ox 


iS x 


2 Jllll^x,2f Xl 4, Jf ig2 f 3 e 4 — JlOll — 


Sx,2JuuS x ,2, 


^hh939i 


— Jiwo 


= Ox 


iS x 


3JllllS x ^a x A, Jf 1 f 2 g 3 e 4 — JllOl = 


S x ,3JllllS x ,3, Jfif 2 f 3 94 — JllW — Ox,iJllHOx,4 



(24) 

So with a state marked, and the four-qubit diffusion transform which is generated by combining two Hadamard 
gates W® 4 with Joooo, a full Grover search for four qubits is available. In principle, if each component of the design 
is available, our scheme would be achievable experimentally. 

Taking the marked state 1/1525364) as an example, a standard quantum circuit for the Grover search al- 
gorithm for N = 4 entries is shown in Fig. 4. The procedure of the Grover search is accomplished with 
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three primitive and sequential steps. First we should perform phase inversion of the desired state by realizing 



J P = 



^/iS2ff3e 4 = Sx,3Sx,2Jii\iSx,2S x ,3- Then we invert all the states with respect to the average ampli- 
tude of all the states, which can be achieved by making use of diffusion transform = Jqqq W® 4 — 

W® 4 a x ^S x fiS x fiS x ,i JuiiSx^Sx^Sx^ax^W® 4 . Finally, a measurement of the whole system is made to obtain the 
marked state 1/1525364). By repeating the first two steps for several times, we can get the marked state with high 
success probability. However, to simplify above considerations, we can reduce the indispensable single-qubit rotations 
(i.e., W, S x , and a x ) to the corresponding transforms Ri and i? i (i = 1, 2 , 3, 4), 



i?i — Ra — 



1 1 

-1 1 



R2 = ^3 — 



-1 1 
1 1 



, R 1 — R 2 — i? 3 — R 4 



1 -1 
1 1 



(25) 



As the qubits 1, 2, and 3 are encoded in the Zeeman sublevels, the transforms Ri and R i (i = 1, 2 , 3) should be 
carried out by a pair of laser beams in Raman process [26], whereas R4 or R 4 operating on the 4th ion could be 
achieved by using only one laser. Along with described in above section, all the necessary operations for Grover 
search are available in a straightforward way. Therefore, based on the reliable control of the laser pulses on stationary 
ions, we could carry out a four-qubit Grover algorithm with four-qubit CPF gates and a series of single-qubit gatings 
according to the operations plotted in Fig. 4. As a result, we consider that our scheme could significantly reduce the 
overhead for QIP tasks with trapped ions due to direct implementation of n-qubit CPF gates. 

As the Grover search involving more than two qubits is carried out only probabilistically, how to obtain an optimal 
search is a problem of much interest. Considering the success probability of the CPF gate (i.e., Eq. (17)) and the 
intrinsic probability of the Grover search itself, we show by numerical simulation the implementation of the Grover 
search in Fig. 5. Due to the similarity, we only demonstrate the search for a marked state in the case of n — 4 and 
5, which demonstrates the decreasing success rates with the increase of m and the qubit number. 



IV. DISCUSSION AND CONCLUSION 



We briefly address the experimental feasibility of our scheme. We may employ S\/2{mj = 1/2), Si/2(rrij = —1/2), 
and D 5 / 2 ( m j — ~ 1/2) of 40 Ca + [27] as the states |/) , \g) and |e) , respectively. The lifetime of D 5 / 2 of the 40 Ca + 
ion is longer than 1 sec, and the transition between |/) and \g) is dipolar forbidden. To carry out the CPF gate 
in the present scheme, the condition 77^™ ax v / 2~7TT„ = it and the approximation relation evf[to / y/2rj] — ► 1 should 
be well satisfied. By choosing suitable values of the (z-axis) trap frequency v and maximal Rabi frequency TO in 
current ion trap experiments, all the above-mentioned relations could be met in the weak-excitation regime Q m <C r\v 
(which implies a negligible AC Stark shift induced by the radiation), where one may expect the unwanted off- resonant 
transition to be disregarded on the precondition that the laser field cannot be made too intense at any rate, which 
thereby is regarded as the dominating limiting factor on the resulting relatively-low gate speed [28]. Taking our 
proposed n-qubit CPF gate as an example, for a typical value r\ = 0.1 of the Lamb-Dicke parameter, we have listed 
the required time to for CPF gate operation in Table II, which clearly indicate that the values of to are prolonged with 
the growth of the number of ions. In addition, the switching rate of the CPF gate will be of order < Z//1000, which, 
from the viewpoint of heating effect, will inevitably bring out disadvantaged impact on the construction of CPF gate. 
Alternatively, in what follows we will take into account another different regimes Ct m -C v/r], where the only new 
requirement we should make is that the laser should take a fixed intensity satisfaying the resonance condition [28] 
&m = v/2, which implies an improvement by two orders of magnitude with respect to the CPF gating time in the case 
of weak-excitation. So to reduce the implementation time to, we may choose such a stronger radiation with Q m = v j2 
for a carrier transition. Our direct calculation based on the model in [28] shows that we could also have Eq. (5) due 
to AC Stark shift and we could should satisfy to ^ 0.574ms, 0.766ms, 0.957ms in the case of n = 3, 4, 5, respectively, 
to achieve our CPF gate. As the heating time of the ground vibrational state of the ions in the linear trap is about 4 
millisec [29], the implementation time in the case of strong radiation seems better in view of avoiding heating. Note 
also that single-qubit operation takes negligible time in comparison with that for many-qubit phase gating, so direct 
calculation shows that one iteration of our proposed Grover search would take 2t ~ 1ms, which is shorter than the 
the heating time of the ions. In this sense, the decoherence originating from the heating effect is not a big obstacle 
for the simulation of Grover search by our proposal when the number of qubit is small. 

TABLE II. List of the required time to for CPF gate and the values of v and Q m (77 = 0.1). 
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regime 


qubit number n 


n = 2 


n = 3 


n = 4 


n = 5 


n = 6 


n = 7 


n = 8 




v (MHz/2n) 


0.250 


0.167 


0.125 


0.100 


0.083 


0.071 


0.063 




n m (kHz/2Tr) 


2.500 


1.670 


1.250 


1.000 


0.830 


0.710 


0.630 




t (ms) 


19.14 


28.71 


38.29 


47.86 


57.43 


67.00 


76.57 




fl m (MHz/2n) 


0.125 


0.084 


0.063 


0.050 


0.042 


0.036 


0.032 




t (ms) 


0.383 


0.574 


0.766 


0.957 


1.148 


1.340 


1.531 



It is generally considered that the computing operation on trapped ions in a linear ion trap would be more and more 
intricate with the increase of the number of ions. Because of the decrease of the spatial separation of the trapped ions, 
individual manipulation is more and more difficult, and meanwhile the vibrational mode spectrum becomes more and 
more unresolved. As a result, the extension of quantum computation from a few qubits to large numbers of qubits is 
quite technically challenging. In this sense, although our proposal is in principle scalable, the currently technical level 
regarding the linear trap restricts the application of our scheme. Nevertheless, eight ultracold ions in the linear trap 
could be individually addressed and entangled [30]. Therefore, if we apply our scheme to these eight ions, a Gorver 
search with 2 8 states could be achieved in a simpler way. 

Alternatively, we may consider the application of our scheme in a multi-zone trap, in which we may carry out our 
scheme on few ions in separate zones, respectively, and then entangle a large numbers of ions by moving the ions 
between different zones [31]. This is a possible way to a large-scale Grover search implementation. 

We have reiterated our scheme to achieve a CPF gate in one step, which could save the implementation time and 
steps compared to conventional methods. Moreover the most impressive feature of our scheme is the higher fidelity and 
higher success rate with more qubits involved, which favors large-scale QIP. Nevertheless, the more qubits, the more 
steps necessary for an optimal Grover search. As the CPF gate proposed here is intrinsically imperfect, a decrease of 
the fidelity is inevitable with more qubits involved, even if we neglect the unpredictable imperfection and decoherence 
in actually experimental situation. Anyway, the curent experimental progress has shown some efficient ways to deal 
with the imperfection [32]. Other basic operations are also available because the qubit encoding in the first n-1 ions 
is like that done in Oxford's group [26] and the qubit encoding in the last ion is employed in Innsbruck's group [27]. 
So we believe that our proposal for Grover search should be available with current or near-future techniques. 

In summary, a potentially practical scheme for performing n-qubit CPF gate as well as the many-qubit Grover 
search algorithm has been proposed in the ion-trap system. We have demonstrated both analytically and numerically 
that, our scheme could be achieved efficiently to find a marked state with high fidelity and high success probability. 
Our proposal could be employed in both the linear ion trap and the multizone trap, and may also be applied to other 
QIP tasks, such as preparation of cluster states for one-way quantum computation [33]. Therefore, we argue that our 
idea in the present paper is not only practical, but also simple and experimentally feasible, which would be helpful 
for working in large-scale QIP devices. 
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Captions of Figures 

FIG. 1. Schematic setup for implementing n-qubit CPF gate and Grover search in a linear trap, where the inset 
shows the ionic level configuration, with bold lines for the states encoding qubits and the arrows for the coupling of 
the lasers to the ions. 

FIG. 2. Infidelity versus the Rabi frequency ratio of the last ion to the others, where the solid, dashed-dotted and 
dashed curves represent the case of n = 3, 6 and 9, respectively. 

FIG. 3. Success probability versus the Rabi frequency ratio of the last ion to the others, where the solid, dashed- 
dotted and dashed curves represent the case of n = 3, 6 and 9, respectively. 

FIG. 4. Quantum circuit of one iteration of the four-qubit Grover search for the marked state \ f1g2g3e4}, where W, 
J^n and 5a; (a x ) are the Hadamard gate, four-qubit controlled phase gate and single-qubit NOT gate, respectively. 
The state of ions is initially prepared in the average superposition state = j(|<7i) + |/i))(|<?2) + |/2))(|<?3) + 

l/3))(lff4) + l e 4))- The operations in the dashed boxes could be reduced to the transforms Ri and R t with i(=l, 2, 3, 
4) denoting the ith ion, which is helpful for experimental implementation. To maximize the search propability, we 
should implement the circuit repeatedly for several times. 

FIG. 5. The search probability for the marked states \f1g2g3e4) and \ f1g2g3e4e5} versus the number of the iterations 
in the case of n = 4 and 5. 
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